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(Dated: May 29, 2007) 

Abstract 

Using a 281 pb _1 data sample collected at the ^(3770) with the CLEO-c detector, we present 
the first absolute branching fraction measurement of the decay D° — > K~ir + ir~ e + u e at a sta- 
tistical significance of about 4.0 standard deviations. We find 10 candidates consistent with the 
decay D° — > K~it + it~ e + v e . The probability that a background fluctuation accounts for this sig- 
nal is less than 4.1 x 10~ 5 . We find B(D° -► R-^^e + v e ) = [2.8^; 4 (stat) ± 0.3(syst)] x 10~ 4 . 
This channel is consistent with being predominantly produced through D° — > iff (1270) e + f e . 
By restricting the invariant mass of the hadronic system to be consistent with If i (1270), we ob- 
tain the product of branching fractions B(D° -► iff (1270) e+z/ e ) • fi(iff (1270) -► A"-7r + 7r - ) = 
[2.5±i;g (stat) ± 0.2 (syst)] x 10^ 4 . Using £(iff (1270) -► if-vr+vr-) = (33 ± 3)%, we obtain 
B(D° -► iff (1270) e+i/ e ) = [7.6±|;J (stat) ± 0.6 (syst) ± 0.7] x 10~ 4 . The last error accounts for 
the uncertainties in the measured iff (1270) — > if~7r + 7r~ branching fraction. 

PACS numbers: 13.20.Fc, 12.38.Qk, 14.40.Lb 
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The understanding of the hadronic mass spectrum in semileptonic decays of charm mesons 
sheds light on non-perturbative strong interaction dynamics in weak decays. In particular, 
an interesting question is whether the charm quark can be considered "heavy," and thus 
theoretical predictions based upon heavy quark effective theory (HQET) can be applied to 
describe some features of its decays. A priori this seems to be an unlikely scenario as, even 
in the Cabibbo-favored transition c — > se + u e , the daughter quark is too light for HQET to 
apply. Nonetheless, this effective theory seems to describe these decays relatively well [1]. 

The decays induced by the quark level process c — > se + v e are dominated by the two 
final states D —>■ Ke + v e and D — ► K*e + v e . CLEO-c has measured exclusive D semileptonic 
branching fractions for all modes observed to date: Ke + i> e , K*e + v e , 7re + z/ e , pe + z/ e , and 
D + — > uje + v e 0, as well as inclusive D — > Xe + v e branching fractions 0]. The sum of the 
exclusive branching fractions and the inclusive branching fractions for D meson semileptonic 
decays are consistent: £B(£>° xcl ) = [6.1 ± 0.2(stat) ± 0.2(syst)]% and £#(£>ici) = [15.1 ± 
0.5(stat)±0.5(syst)]% while B{D° -> Xe + u e ) = [6.46±0.17(stat)±0.13(syst)]% and B{D+ -> 
Xe + v e ) = [16.13 ±0.20(stat) ±0.33(syst)]%. Nonetheless, there is some room left for higher 
multiplicity modes. 

The quark model developed by Isgur, Scora, Grinstein, and Wise Q, later updated to 
include constraints from heavy quark symmetry, hyperfine distortions of wave functions, 
and form factors with more realistic high recoil behavior [lj], is the only one to provide 
quantitative predictions for the partial width of decays such as D — > i^i(1270) e + v e . In 
general, we expect the decay mediated by the quark level process c — > se + v e to be dominated 
by the ground state pseudoscalar and vector daughter mesons. The low available phase space 
makes it less likely to produce heavier mesons, such as P-wave or first radial excitations of 
the su and sd quark states. The lightest excited state is the 1^(1270). This model predicts 
that the partial width T(D — > f^i(1270)e + z/ e ) is 2% of the total T(c — > se + v e ), and that 
decays to other excited resonances are suppressed by at least a factor of 10 more. 

Little is known about D° — ► i^i(1270)e + z/ e to date. The fixed target experiment E653 
Q reported a 90% confidence upper limit of B(D° -> R-ix+tj- p, + v^) < 0.037 x B(D° -> 
K~/j i + i/ ix ). This Letter is the first report on a signal for the decay D° — ► K~7T + 7T~e + i> e . 
We use a 281 pb _1 data sample collected at the ^(3770) with the CLEO-c detector 
0]. The three major subsystems of this detector are the charged particle tracking cham- 
bers, the Csl electromagnetic calorimeter, and a Ring Imaging Cherenkov (RICH) charged 
particle identification system. All these components are critical to an efficient and highly 
selective electron and positron identification algorithm. The Csl calorimeter measures the 
electron and photon energies with an r.m.s. resolution of 2.2% at E = 1 GeV and 5% at 
i?=100 MeV. One of the key variables for e identification, E/p, uses E, the energy measured 
in the calorimeter and p, the momentum measured in the charged particle tracking system. 
The tracking system is composed of a 6-layer inner drift chamber and a 47-layer main drift 
chamber. The main drift chamber also provides specific ionization (dE/dx) measurements 
for charged particle identification. In addition, charged particles are identified using the 
RICH detector [§[. Combining information from these detector subsystems, we achieve effi- 
cient and selective charged particle identification over the entire momentum region relevant 
for the decays studied. 

We use a tagging technique similar to the one pioneered by the Mark III collaboration 
Details on the tagging selection procedure are given in Ref. [10|]. We select events containing 
a fully reconstructed D° -> K + n~, D° -> K+tv-jt , D° -> K+n^n+n-, or D° -> i^TT+TT 
decay, which we call a tag. (Mention of a specific mode implies the use of the charge 
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conjugate mode as well throughout this Letter). Two kinematic variables, namely energy 



difference, AE = E tSug — -Ebeam, and beam-constrained mass, Mb c = yE^ eSbm /c A — |pta g | 2 /c 2 , 
are used to select tag candidates, where -Ebeam represents the beam energy and (E tag ,p tag ) 
represent the 4-vectors of the D° tag candidate. We first require \AE\ to be less than 0.020 
to 0.030 GeV, depending upon the mode considered. Figure [1] shows the Mb c spectra for 
events that satisfy the \AE\ requirement for the four tagging modes considered. In order 
to determine the total number of tags, we fit the M^ c distribution with a signal shape 



composed of a Crystal Ball function [ll| and a Gaussian, and an ARGUS function [1J 



parameterizing the background in the fit. The signal window is chosen as 1.858 GeV/c 2 < 
Mbc < 1.874 GeV/c 2 . In order to extract the tag yield, we integrate the signal shape 
within this M^ c interval. Alternatively, we count tag candidates in the Mb c signal window 
and subtract the combinatorial background obtained by integrating the background function 
from the fit. The total number of tags obtained with the former method is [257.4±0.7(stat)] x 
10 3 ; the second method gives [257.7 ± 0.6(stat)] x 10 3 . The agreement is excellent and we 
use the latter number as the total number of tags in our sample. The difference between 
the two tag yields is included in a systematic uncertainty. 

In each event where a tag is found, we search for a set of tracks recoiling against the tag 
that are consistent with a semileptonic decay. We select tracks that are well-measured and 
have a helical trajectory approaching the event origin within a distance of 5 cm (5 mm) along 
the beam axis (in the plane perpendicular to the beam axis). Each track must include at least 
50% of the main drift chamber wire hits expected for its momentum and have momentum 
greater than 50 MeV/c. We search for a positron among well reconstructed tracks having a 
momentum of at least 200 MeV/c, as the electron identification becomes increasingly difficult 
at low momenta. We also require | cos#| < 0.90, where 6 is the angle between the positron 
direction and the beam axis. The positron selection criteria are discussed in Ref. 0]. They 
have an average efficiency of 95% in the momentum region [0.3 — 1.0] GeV/c, and 71% in the 
region [0.2 — 0.3] GeV/c. In addition, we search for a good track consistent with a K~ and 
two oppositely charged tracks consistent with pions. Hadron track identification criteria rely 
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on dE/dx information from the drift chamber for tracks with p < 0.7 GeV/c. For tracks 
with p > 0.7 GeV/c, in addition to dE/dx measurements, information from the RICH 
detector [8j] is used to improve the K-tx discrimination. In the momentum range relevant for 
this analysis the K-tt misidentification probability is negligible. The e-n misidentification 
probability, determined experimentally with radiative Bhabhas, has an average value of 17% 
for electron momenta below 0.2 GeV/c, and is about 1% for higher momenta. 

As the decay mode that we are investigating is rare, efficient background suppression 
is critical to achieve adequate sensitivity. Accordingly, we require that only four charged 
tracks be present in the event in addition to those used in the tag reconstruction. The 
dominant source of background in this analysis arises from events in which the detected 
positron comes from a 7 conversion (7 — > e + e~), or a tt° Dalitz decay (7r° — ► e + e~ r y). 
This background is equally likely to produce (e + K~) combinations, which we call right-sign 
events (RS), and {e~K~) combinations, which we call wrong-sign events (WS). Typically, an 
e + e~ pair arising from a conversion 7 or a tt° Dalitz decay has a strong angular correlation 
with almost collinear angular orientation of the two particles. For signal events, the opening 
angle between the e + 7r~ pair tends to be large. We therefore include a requirement that 
the opening angle be greater than 20°. This requirement eliminates most of the background 
from conversion 7's or 7r° Dalitz decays, while reducing the signal efficiency by only 1.7%. 

In this semileptonic sample, signal candidate events are selected using the missing mass 
squared MM 2 defined as 

MM 2 = (£ beam - £ Etf/c 4 - (-p tag - Y,Pif/c\ (1) 

i=l i=l 

where p ts , g is the momentum of the fully reconstructed tag, and (Ei,pi) represent the en- 
ergy and momentum of the four tracks in the D° candidate. For signal events the MM 2 
distribution is centered at zero, as it represents the invariant mass squared of the missing 
u e . According to Monte Carlo simulation of our signal semileptonic channel, the MM 2 dis- 
tribution has a resolution consistent among the tag modes with a standard deviation (a) 
of 0.00594 ± 0.00010 (GeV/c 2 ) 2 . Figure M shows the measured MM 2 distribution for RS 
events in the data as well as the estimated background, derived from GEANT-based Monte 
Carlo simulation 13[ in combination with particle misidentification probabilities derived 



from data. In addition we estimate the background directly from the WS events in data. 
We define a signal window as |MM 2 | < 0.02 (GeV/c 2 ) 2 . There are 10 events in the signal 
window of MM 2 as shown in Figure [2J 

Another interesting observable is the invariant mass of the K~tt + tt~ hadron system. 
Figure [3] shows the invariant mass of the K~7r + 7r~ system for RS candidate events, compared 
with the expectation from the ISGW2 model [l[, which provides the best representation 
of our data, where the hadronic system forms the Xi(1270) resonance. The measured 
distribution is in reasonable agreement with this model. 

We have performed several studies to determine possible background sources. A Monte 
Carlo sample incorporating all the information available on D meson decays and 40 times 
bigger than our collected data demonstrates that the dominant background comes from 
conversion 7's or 7r° Dalitz decays. As the e to tt misidentification probability may not 
be modeled accurately by our Monte Carlo simulation, the background from Dalitz decays 
is evaluated by folding the e spectra from simulated D° — > K^7i + tt° decays with the e to 
it misidentification probability derived from a radiative Bhabha data sample. This study 
predicts that 1.56 ± 0.22 background events are due to this source if no requirement on the 
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FIG. 2: Missing mass squared (MM 2 ) distribution for the RS sample D° — > K~i: + T:~e + v e . The 
dashed histogram represents the estimated background. Events with MM 2 within the two arrows 
are considered signal candidates. 
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FIG. 3: Invariant mass of the hadronic system in the data for D° — ► K~7t + Ti~ e + v e . The dashed 
histogram represents the predicted distribution obtained using a Monte Carlo simulation according 
to the ISGW2 model, assuming all the K~it + it~ are If i (1270) decay products. The region within 
the two arrows defines the invariant mass range used to select the iff (1270) resonance. 



K 7r + 7r invariant mass is applied. A study of the WS data gives one background event, in 
agreement with the previous estimate. In addition, there are small background components 
from the decays D° -> R-n+n+n- (0.2 ± 0.1) and D° -> K-n + n + n-n° (0.1 ± 0.1), both 
estimated with Monte Carlo samples. We have also studied non-DD contributions at this 
center-of-mass energy, such as those from the continuum (e + e~ — > qq, where q is a u, d, 
or s quark), radiative return production of ip(2S), and e + e~ — > r + r _ processes, and we do 
not find any background from these sources. Summing up all background contributions, we 
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find that 1.86 ± 0.25 (stat) events are consistent with background. We have also studied 
this sample with a requirement on the invariant mass of the Kitit system optimized for 
the decay D° — > K{ (1270)e + z/ e , using the variable S/y/S + B, where S is the number of 
signal events predicted from Monte Carlo simulations and B is the number of estimated 
background events. For the optimal invariant mass interval, [1150 — 1500] MeV/c 2 , we find 
8 candidate events and an estimated background of l.Oifj^stat) events, with no events in 
the WS sample. 

The reconstruction efficiency depends on the invariant mass of the fT~7r + 7r~ system 
(Mhad). A larger fraction of the electron spectrum is below the momentum cut of 0.2 
GeV/c for higher Mhad, and the spin and parity of the final hadronic state influence the 
electron spectrum shape as well. For example, the ISGW2 model studies all the P-wave su 
and du hadronic final states, as well as the corresponding radial excitations. Among the 
P-wave states, the 3 ^ 2 Pi are identified with the i^i(1270), and 1 ^ 2 P 1 are identified with the 
i^i(1400). The latter has a much softer electron spectrum, and therefore our efficiency for 
detecting it is smaller. We have studied the signal reconstruction efficiency with the ISGW2 
model, including different mixing percentages of the 3 / 2 Px and l / 2 Pi final states, as well as 
a phase space model for the distribution of the Mh a d- With the Monte Carlo simulation 
based on the ISGW2 model, we obtain e = (10.78 ± 0.23)% for the full Mh a d range and 
e = (10.53 ± 0.22)% with the 1^(1270) mass requirement ([1150 - 1500] MeV/c 2 ). 

The absolute branching fraction for D° — > K~7T + 7T~e + b> e is obtained using B = (iV s — 
-WbV(eeff-Wtag); where N s is the number of signal events, N b is the number of background 
events, N taig is the number of tags, and e e fj is the effective efficiency for detecting the semilep- 
tonic decay in an event with an identified tag. This effective efficiency includes a correction 
term C = e^ g /e tag accounting for the small difference in tag reconstruction efficiency in 
events containing the semileptonic signal and in generic DD events. The average value of 
C is 1.036. We obtain B(D° -> R-tt+tx- e + u e ) = (2.8±\i ± 0.3) x 10~ 4 , without applying 
any invariant mass requirement. If we apply the ^(1270) invariant mass requirement, we 
obtain B(D° -> K{ (1270)e+z/ e ) • £(7^(1270) -> R-n+n-) = (2.5±i;g ± 0.2) x 10" 4 . The 
smaller systematic uncertainty is derived by the fact that the model dependence can simply 
be estimated by varying the form factors in the ISGW2 model. In this case we do not 
need to model a broader invariant mass distribution for the K~7t + 7t~ system. Note that the 
probability for 1.86 background events to fluctuate to 10 or more events, taking into account 
a 0.25 event Gaussian uncertainty, is 4.1 x 10~ 5 , corresponding to a significance of about 4.0 
er. The result with the K^-k + -k~ mass requirement has similar statistical significance. 

The systematic uncertainties for the branching fractions are listed in Table [I] and are 
quoted as relative to the measured branching fraction. The uncertainty on the tag yield is 
estimated from varying the background functions. Systematic uncertainties on track finding 
and hadron particle identification efficiencies are reported in Ref. while electron iden- 
tification efficiency is reported in Ref. [jjj]. The sensitivity to the requirement on the e + 7r~ 
opening angle has been evaluated by repeating the analysis after changing the requirement 
by ±5°. The model dependence of the efficiency is studied using an alternative invariant 
mass distribution for the hadronic system governed by phase space. In the analysis where we 
apply a mass requirement on the K~-k + -k~ system, the model dependence of the efficiency 
is estimated by varying the form factors in the ISGW2 model, and the corresponding uncer- 
tainty is found to be 4%. The background uncertainty is derived by changing the measured 
fake probabilities within their errors. 

In summary, we have presented the first measurement of the absolute branching fraction 
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Systematic errors (%) 



Number of tags 
Tracking 

PID Efficiency (hadrons) 

PID Efficiency (electrons) 

Opening angle cut 

M Kw7T cut 

Model dependence 

Background 



10.0 
5.3 



0.5 
1.3 
1.9 
1.0 
1.5 



0.5 
1.3 
1.9 
1.0 
1.5 
1.7 
4.0 
5.3 



Total 



11.9 



7.5 



TABLE I: Systematic errors on 7r + 7r e + v e branching fraction. The first column applies 

to the analysis without iG(1270) mass cut, the second to the analysis with the K \ (1270) mass cut 
as described in the text. 

B(D° -> K-7T + 7v-e + u e ) = [2.8±\i (stat) ± 0.3 (syst)] x 1(T 4 . The invariant mass of the 
hadronic system recoiling against the e + v e pair is consistent with iff (1270). By requiring 
A^had to be within the [1150 — 1500] MeV/c 2 mass window, we obtain the product branching 
fraction B(D° -> K{ (1270)e + z/ e ) • (£(A7(1270) -> K~n + n~) = [2.5±i;g stat ± 0.2] x 10" 4 . 
The statistical significance is about 4.0 standard deviations. Using the iCf (1270) decay 
modes reported in the PDG [13], we calculate the iff (1270) — > K~7t + tt~ branching frac- 
tion to be (33 ± 3)%. Thus, the absolute branching fraction B(D° — > K{ (1270)e + i/ e ) is 
[7.6lJ;o (stat) ± 0.6 (syst) ± 0.7] x 10~ 4 . The last error accounts for the uncertainties in 
the measured ITf (1270) branching fractions. This channel is found to be 1.2% of the to- 
tal semileptonic width. The ISGW [H model predicts this fraction to be about 1%, while 
the ISGW2 model [lj predicts this fraction to be about 2%; hence the measured branching 
fraction and K~7t + tt~ invariant mass are consistent with quark model calculations. 
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